RADEMACHER'S INFINITE PARTIAL FRACTION 

CONJECTURE IS (ataost certainly) FALSE 

ANDREW V. SILLS AND DORON ZEILBERGER 

"It is a capital mistake to theorise before one has data. Insensibly one begins to twist facts to 
suit theories, instead of theories to suit facts." — Sherlock Holmes to Dr. Watson [2l p. 63]. 

Abstract. In his book Topics in Analytic Number Theory, Hans Rademacher 
conjectured that the Umits of certain sequences of coefficients that arise in the 
ordinary partial fraction decomposition of the generating function for parti- 
tions of integers into at most A'^ parts exist and equal particular values that 
he specified. Despite being open for nearly four decades, little progress has 
been made toward proving or disproving the conjecture, perhaps in part due 
to the difficulty in actually computing the coefficients in question. In this pa- 
per, we provide a fast algorithm for calculating the Rademacher coefficients, 
a large amount of data, direct formulas for certain collections of Rademacher 
coefficients, and overwhelming evidence against the truth of the conjecture. 
While the limits of the sequences of Rademacher coefficients do not exist (the 
sequences oscillate and attain arbitrarily large positive and negative values), 
the sequences do get very close to Rademacher's conjectured limits for certain 
(predictable) indices in the sequences. 



Iimportant Note 

This article is accompanied by the Maple package HANS, downloadable from 
http : //www .math . rutgers . edu/~zeilberg/tokhniot/HANS . 

The "front" of this article, 
http : //www. math. rutgers . edu/'zeilberg/mamarim/mamarimhtml/hans .html 

contains lots of supporting input and output files. 



1. Introduction 

Let pAr(n) denote the number of partitions of the integer n into at most N parts 
The generating function of pAr(n), 



^ 1 



n>0 j=l 
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may be decomposed into partial fractions: 

N N [N/k] 

(11) TT = V V V ^''"''^^ 

j=l k = l 0<h<k 1=1 ^ ' 

gcd{h,k) = l 



We shall refer to the Ch.k,i{N) defined by (1.1) as the Rademacher coefficients. 

Near the end of his posthumously published masterpiece Topics in Analytic Num- 
ber Theory p. 302], Hans Rademacher made the following conjecture; 

Rademacher's Conjecture. For all integers h, k, I such that < h < k, gcd{h, k) ■ 
1 and I > 1, ImiN^oo Ch,k,iiN) exists and equals 

(1-2) R.^k, := -2n {-) ^3/2 (-^(« + 

evaluated at a ^ ^, where s{h, k) = J2l=\ - J ' ^) {'f - l'f\ " 5) ihe 
Dedekind sum, is the forward difference operator, so that 



A^„/(a) = £(-l)'^(^)/(a + j-/.), 

h=0 



una 

^3/2(V)--^(2cos(2,)-^). 

If the Rademacher conjecture would have been true then it would have followed 
that 







- 


1 + 


25 \^ 


144 




1225 


+ 4 




57r 


1 cos 






I2 




= -0.292927573960...), 



lim Co 1 2(^) = -Ro 1 2(= + = 0.1897670688440 



N 

lim Ci,2.i(A^) = i?i.2.i(=-^ fcos^-^sin^) = 0.093882853484...). 
-/v-i-oo 25 \ 12 57r 12/ 

Remark 1.1. The floating-point approximation for the value of i?o,i,i stated by 
Rademacher [H p. 302] was erroneous, as were the exact values of iio,i,2 and i?i,2,i 
(for the latter he gave exactly one half of the correct value). These erroneous values 
were quoted, without correction, by Andrews W, p. 



Rademacher supplied (with one error) a table of values for Co,i,i{N), Co,i,2{N), 
and Ci^2,i{N) for = 1, 2, 3, 4, 5, and in fact these values are not too far off from 
his conjectured "A — 00" cases. 

Rademacher began work on the book in which this conjecture appeared [5] no 
later than 1944 and was still working on it at the inception of his final illness. 
Thus as the final version was edited and published by Rademacher's students Emil 
Grosswald, Joseph Lehner, and Morris Newman, after Rademacher's death in 1969, 
we will never know whether Rademacher came to doubt the truth of the conjecture 
after he had written it down. However, George Andrews reports that Rademacher 
discussed the conjecture in a course he taught at the University of Pennsylvania 
during the 1961-1962 academic year. 
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In [4], Augustine Munagi considered a different type of partial fraction decom- 
position called g-partial fractions, and proved a special case of the analog of the 
Rademacher conjecture, relative to the q-partial fraction decomposition. 

We should also note that the first to cast doubts on the Rademacher conjecture 
were Jane Friedman and Leon Ehrenpreis. Ehrenpreis [3, p. 317] stated, "If one 
attempts to carry out the usual type of partial fraction decomposition of the parti- 
tion function term-by-term, it is difhcult to compute the coefficients. My student, 
Jane Friedman, spent a great deal of time trying to apply computer algorithm 
methods to compare the coefficients with those of Rademacher. . . Unfortunately, 
the computer study proved inconclusive." 

In this article, we present overwhelming evidence against this conjecture, taken 
literally, but we will present ample evidence for a modified conjecture. Additionally, 
we will present a fast algorithm for generating the Rademacher coefficients, and 
formulas for a selection of particular Rademacher coefficients. 

2. Empirical evidence against the Rademacher conjecture 

2.1. The actual behavior of the sequences Co,i./(iV), I — 1,2,3, ... . At 
http : //www .math . rutgers . edu/~zeilberg/mamarim/mcmiarimhtml/hans . html , 

there are links to various files including 

• the sequences Co,i,z(A^) for 1 < / < 10 and 1 < iV < 850, in Maple-readable 
format, 

• the sequences of floating point approximations to Co,i,;(Af) for 1 < ^ < 40 
and 1 < iV < 1000, in Maple-readable format. 

Figures IT] and [2] are graphical summaries of Co,i,i(iV). 
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Figure 1 


Graph of Co,i,i(iV) for N from 1 to 100, together with 


the line y 


— Ro.lA 



Note that Coa,i(25) differs from i?o,ia, Rademacher's conjectured value of Co, 1,1(00), 
by less than 0.000032, but then things go down hill (for the conjecture) from there, 
and get particularly bad after about n — 150. Numerical evidence points to the 
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Figure 2. Graph of Co,i,i(^) for N from 1 to 200. 



sequence Co,i,i(-/V) oscillating and attaining arbitrarily large positive and negative 
values. The same thing is true for other sequences Ch,k.i{-^)'^ see the output file at 
http : //www. math. rutgers . edu/~zeilberg/tokhniot/oHANSll . 

By examining the graphs in Figures [l] and [2] and the associated numerical data, 
it seems reasonable to state the following alternative conjecture: 

Conjecture 2.1. Co,i,i(A^) is an oscillating function of N of "period" 32, with local 
maxima (resp. local minima) that attain arbitrarily large positive (resp. negative) 
values as N increases. 

By "period 32" we mean that the peaks and valleys, eventually, recur at a period 
of 32. We also noticed, numerically, that the elevations and depths of successive 
peaks and valleys roughly grows exponentially with a factor around 8. Specifically, 
Co,i,i (A^) has local maxima at iV = 3, 4, 33, 66, 99, 131, 163, 195, 227, 259, 291, 323, 
355, 387, 419, 451, 483, 515, 547, 579, 611, 643, 675, 707, 739, 771, .... The ratio of con- 
secutive local maxima is 

{1, 1.103504574, 0.6965131681, -0.7709983810, 13.63072659, 6.485614677, 
6.289519948, 6.547018652, 6.785098547, 6.992410281, 7.161220864, 7.301859590, 
7.420337150, 7.521483398, 7.608822684, 7.684977203, 7.751953124, 
7.811301903, 7.864245038, 7.911756412, 7.954622120, 7.993483579, 

8.028869316, 8.061218737, 8.090900135, . . . } 

C'o,i,i(-^) has a local minimum half way between each local maximum, at = 
18, 50, 83, 115, 147, 179, 211, 243, 275, 307, 339, 371, 403, 435, 467, 499, 531, 563, 595, 
627,659,691,723,755,787,.... 

More precise conjectured asymptotics for Co.i_;(A^) for I between 1 and 15 can 
be gotten from 

http : //www .math . rutgers . edu/~zeilberg/tokhniot/oHANS10 
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It appears that all Co,i,/(A^) have a "period" of 32. The locations of the local 
maxima and minima of Co,i,;(iV), from = 99 until N — 803 occur when = 3 
(mod 32) and ee 19 (mod 32) respectively. For N > 803, they become N = 2 
(mod 32) and iV = 18 (mod 32) respectively. This gives (very meager) evidence of 
a "shifting of the perihelion" , but one would need to go much further to investigate 
this. More generally, for Co,i,;(A^), and for 100 < N < 800 these locations are 
congruent to 12 — 9^ (mod 32), and 28 — 91 (mod 32) respectively. Probably these 
too would eventually get shifted (ever so slowly). 

2.2. How to compute the sequences Co,i^i{N) fast. If you use the definition of 
^^0,1,1(^)1 or Andrews's formula [TJ p. 388, Theorem 1], you can't go very far, even 
with Maple. Rademacher calculated Co,i,i(iV) for N = 1,2,3,4,5, presumably by 
hand, and made an error in the = 5 case. Andrews [1] p. 388], who had access 
to a computer algebra system in 2003, corrected Rademacher's error at A^ = 5 and 
extended the list to A^ = 6, 7, 8. 

We need to be more clever. A fast recurrence for Co,i,;(A^) can be derived as 
follows. Since 

fr 1 ^^ Co.i.z(iV) , 
we can multiply both sides hy [x — 1)^ and get 

N N-1 

{x-^r\lY^-Y.Dr{N){x~ir^... . 

j = l r=0 

Once we know Dr{N), we can find Co,i,i(Af), since they equal D]si-i{N). It remains 
to find a fast recurrence for Dr{N). 

By definition, we have: 



1 — X 



^ _ ^ DriN)ix -1Y\=}^ Dr{N ~ l){x - ly 

\r=0 / r=0 

Letting z — x ~ 1^ this is: 

\r=0 / r=0 

By the binomial theorem, 

\a=0 ^ ^ / \r=a J r=0 

Equating coefficients of z"^ we get: 

^ / A^ 



And finally: 



a=l 



NDr{N) + ^[^^]Dr-a{N)=-Dr{N-l) 



a— 1 ^ 

This is implemented in procedure C01(1,N) of HANS. 



6 



A. V. SILLS AND D. ZEILBERGER 



The same argument leads to efficient recurrences for Ch.k,i{N), except tliat 
now we liave to distinguislr between tlie case wlren N is divisible by k and when 
it is not, yielding two different recurrences. This is implemented in procedure 
ChklN(h,k,l,N) of HANS. 



2.3. Ci^2AN). The values of Ci,2a(A^) for 1 < < 700, in both exact rational 
form and approximate floating point form are provided at 

http : //www . math . rutgers . edu/~zeilberg/tokhniot/oHANS3 . 

Graphical summaries are provided in Figures [3] and |4] As the graphs make 
clear, it is better to regard Ci^2.i[N) as two separate subsequences, Ci_2,i(2n) and 
(^1,2,1 (2ri + 1). For each subsequence, we see behavior that is similar to that of 
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Figure 3. Graph of Ci,2,i(A^) for iV from 1 to 150, together with 
the line y — i?i,2,i 



3. "Top down" formulas for the Rademacher Coefficients 

3.1. Co,i,/(A^). As Rademacher already pointed out, it seems hopeless to get a 
closed-form formula for Ch,k,i{N) for ^ = 1, 2, . . . , hut if you work your way down 
from the "top" , one can conjecture, and then rigorously prove explicit formulas for 
Ch,k,N-r{N)^ that alas, get increasingly more complicated as r gets larger. 

Conjecture 3.1. 

where, for r > 0, Po.i.TV-rC-^) *s « convex, alternating, monic polynomial of degree 
2r whose only real roots are and 1. 



Theorem 3.2. Explicit formulas for the Po,i^N~riN) of Conjecture 3.1 may he 
given for any specific r, in particular, we have 

(3.1) Po,i,NiN) = 1 
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Figure 4. Graph of Ci.2,i(iV) for N from 1 to 300, together with 
the Une y = i?i,2,i 



(3.2) 



Po,i,N-iiN) = - N 



(3.3) 



Po.i,N-2{N) = N'' ■ 



22N^ 13N^ 26N 



/,rN .rfi 13^^ 437V4 ^^^^^ 98iV2 56A^ 
(3.4) Po,i,N-3{N) = ^ + ^^ 25iV3 + — — 



(3.5) Po,i,N-4{N) ^ N"" 



20N'^ 862iV6 21104iV5 29039Af'' 14548A^3 



27 



225 



135 45 
9892iV2 42896iV 



27 



225 



Remark 3.3. Readers desiring formulas for Co,i,Ar-, (Af) for r > 4 are directed to 
the ChkFormula procedure in the HANS Maple package. 

Remark 3.4. From the above (and additional data not reproduced here but available 
at the website), we may deduce that 



200r6 - 300r5 + 40706r4 + 42939r3 - 257509r2 + 173964r 2r-3 



150 • 93 



lower degree terms. 



For s > 1, the coefRcients of N'^'^~'' again appear to be polynomials in r of degree 
2s, whose real roots include and 1, although they are not convex and may have 
additional real roots. 
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Proof of Theorem \3.S\ Define Gjv := {x — 1)^ Fm{x). Then Gat has a Taylor series 
expansion about x = 1, whose first N coefficients are the Rademacher coefficients: 



N-l 



Clearly, 
(3.6) 



Gn = E ^o,i,Af-j(^)(2: — 1)-' + higher degree terms . 

3=0 



{1-x'')Gn^{x-1)Gn-i. 



Expanding (1 — x^) on the left hand side of (3.6 1 as a Taylor polynomial about 
X = 1, we have 



(3.7) 



N 

E 



N 



' N-l 



ly I ^ Ca^i,N-j{N){x — ly + higher degree terms 

E ^o.i,N-j{N — l)(a; — l)"* + higher degree terms 



Comparing the coefficients of {x — 1)^ on both sides of (3.7), we find 
-iVCo,i^Ar(iV) = Coa,jv-i(iV-l). 
Solving the recurrence with the initial condition Go,i,i(l) — —1, yields 



(3.8) 



which is (3.1 1. 



Comparing the coefficients of [x— 1)^ on both sides of (3.7), we find, taking into 

N\ (-1)^ 



account (3.8 1, 
(3.9) 



l.JV-l 



[N) 



Co.l,N-2{N -I) 



with initial condition (70,1,1(2) — — j yields 



(3.10) 



(-1)^+1 
A{N -2)V 



which is (3.2 ). 

Comparing the coefficients of [x— 1)^ on both sides of (3.7), we find, taking into 
account (3.8| and (3.10), 

fN\ (-1)^+1 fN\ (-1)^ 
(3.11) ~ A.Go,i,._.(iV) - ( J 1^ - ( 3 ) SvT = - 1) 



with initial condition (^0,1,1(3) — —^2 yields 



(3.12) 



C'o,l,Ar-2(^) 



(-l)^(9iV2 - 13A^ + 26) 
288(iV-2)! 



which is (3.3). Results for larger r follow analogously. 



□ 
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3.2. Ci^2,i{N). Let us now define Gat analogously to Gat. Let 

Then Gjv has a Taylor series expansion about x = —1, whose first [^J coefRcients 
are the Rademacher coefficients: 

LW/2J-1 

Gn — ^ Gi,2,L7v/2j-r(^)(2; + 1)'^ + higher degree terms . 

We now abandon the use of the floor function. Notice that 

(3.13) (1 - x2"-i)(l - a;2")G2„ = (x + l)G2„-2. 

Thus, by expanding the two left most factors on the left side as a Taylor series 
about X — ~1, 



^ ^471- 1\ /2n\ (In - 1 



T ) \ 



r \ r 



{x + 1)'^ 



(3.14) <j ^ {-ir+^ 

^ Gi,2,n-r(2n)(x + 1)*" + higher degree terms ] 

n-l 

— ^ Gi,2,n-r(2n — 2){x + ly + higher degree terms 



r=l 



By comparing coefficients of (a; + 1)'' in both sides of (3.14) and solving the re- 
currences, we obtain formulas for Gi^2,n-r(2Ti) analogous to those for Go,i.Ar-r(-^)- 

(3.15) Gi^2,„(2n) = ^ 



22"n! 

71 

(3.16) Gi,2,„-i(2n) = 



(3.17) Gi,2.„-2(2n) 



22"(n- 1)!' 
18^3 - + 15n + 2 



9-22«+2(„-l)! ■ 

Of course, the observation (1 — a;2")(l — 2:2"+^)G2ri+i = {x + l)G2n-i leads to 
analogous formulas for the Gi,2,ii-r(2n + 1), e.g., 

(3.18) Gi,2,„(2n + 1) ^ 

(3.19) Gi,2,„-i(2n+l) = 

(3.20) Gi,2,„-2(2n+l) 



22«+i „!' 

2^2 + 2n + 1 
22"+2 n! ' 

ISn^ + 46n'* + 61n^ + 53^2 + 29n + 9 



9-22»+3 (n + 1)! 

Clearly, the same idea can be used to find formulas for Ch.k,n-j{kn + r) for 
any /i, fc, j', r. This has been implemented in the procedure ChkFormula in the 
HANS Maple package. For those desiring automatically generated papers, containing 
both formulas of this type and their proofs, please use the HaoisTopDownAutoPaper 
procedure in the HANS Maple package. 
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4. Close Encounters of the Rademacher Kind 

While it appears that limjv->.oo Ch,k,i{N) does not exist for any {h,k,l), we can 
nonetheless define Bh^k,i to be the N which comes closest to Rh,k,l- This is imple- 
mented in the CloseEncounters procedure in HANS. 



I 




C'o,i,((-Bo,i,() — -^0,1,' 


Co,i,((-So,i,;)/-Ro,i,( 


1 


25 


0.0003177 


0.99989 


2 


47 


0.0001434 


0.99924 


3 


71 


0.0000828 


0.99991 


4 


149 


0.0000009 


1.00001 



Notice that the first few values of 5o,i,; are close to 24Z. This motivates us to 
consider comparing Co^i^;(24^) to i?o,i.i- 



I 


|Co,i,;(240-i?o,i,i| 


|Co,i,j(240/i?o,i,j| 


1 





0053741095 


1.018346206 


2 





0015044594 


1.007927400 


3 





00033240887 


0.996241370 


4 





00004427030 


1.001376635 


5 





000011288321 


0.9988220859 


6 





000001686611 


1.0006971253 


7 





0000001275687 


0.9997575030 


8 





0000000110523 


1.0000986383 


9 





00000000239242 


0.9999562770 


10 





000000005333208 


1.0000141594 


11 





0000000187490584 


0.9999947242 


12 





0000000393434274 


1.0000017401 



Thus we have some evidence that even though the A'' for which Co,i,i(A?') is 
closest to Ra.i.i is not = 24Z, Co.i,((24Z) seems to provide a good approximation 
to -Ro,i,i) a-nd the approximation seems to be improving as I increases. 
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